Die kontroUierte Kettenreaktion 
(The controlled chain reaction) 

Andreas Aste 

Paul Scherrer Institute, CH-5232 Villigen PSI, Switzerland 
February 24, 2012 



On 



Abstract 

Es werden die Grandlagen der kontroUierten nuklearen Kettenreaktion in kompakter Form diskutiert. 

<N 

Einfiihrung 

<N 

^ Der Inhalt der vorUegenden Schrift befasst sich mit der quantitativen Beschreibung der Kettenreaktion, wie 
a sie in (kommerziellen) Leistungsreaktoren oder Forschungsreaktoren zu steuern ist. Das eigentliche Prinzip 
^ der Kettenreaktion ist simpel: die Spaltung von ^'^^U-, ^'^^Pu- oder ^'^'^U-Kernen durch moderierte, d.h. nieder- 
energetische Neutronen setzt wiederum Neutronen frei, welche potenziell weitere Spaltungen im Brennstoff des 
Reaktors auslosen. Im Gegensatz zur unkontroUierten Kettenreaktion, bei welcher hochenergetische Neutronen 
I (— ( ' im komprimierten, hochangereicherten Spaltstoff einer Kemwaffe innert kiirzester Zeit (10~^s) einen grossen 
Anteil der vorhandenen Kerne spalten und deren Spaltenergie freisetzen, ist die Kettenreaktion im Reaktor 
dank dem Vorhandensein verzogerter Neutronen steuerbar. Diese Neutronen werden nicht unmittelbar bei der 
^ Kemspaltung freigesetzt, sie werden vielmehr mit einer durchschnittlichen Verzogerung von etwa 13 Sekun- 
, den durch die bei der Spaltung entstehenden Kembruchstilcke emittiert. Ebenso wie die bei der Kernspaltung 

prompt freigesetzten Neutronen miissen auch die verzogerten Neutronen durch ein Moderatormaterial abge- 
bremst werden, damit sie mit ausreichender Wahrscheinlichkeit weitere Kemspaltungen auslosen konnen - 
schnelle Neutronen mit Geschwindigkeiten im Bereich von einigen tausend Kilometem pro Sekunde spalten 
^ die Spaltkeme nicht sehr effektiv. 

in 

^ Moderierte Neutronen bewegen sich in einem Reaktor bei Temperaturen um 300 C mit einer charakteristis- 

Q chen Geschwindigkeit von etwa 3km/s. Auch wenn diese Neutronen fortwahrend KoUisionen mit Atomkemen 

CN erleiden und daher jeweils nur geringe Distanzen der Grossenordnung eines Zentimeters von Stoss zu Stoss 

^ zuriicklegen konnen, so breitet sich eine lokale Storung der Neutronendichte innerhalb eines Reaktors von 

^ wenigen Metem raumlicher Ausdehnung innerhalb Bruchteilen einer Sekunde aus. Dieses Verhalten rechtfer- 

^ tigt weitestgehend eine koUektive Beschreibung der Neutronen durch das sogenannte Punktmodell. Ausgehend 

H von einem Minimum an Annahmen und Daten konnen im Rahmen dieses Modells zahlreiche interessante Aus- 
sagen iiber die Kinetik der Neutronenbilanz gemacht werden. 
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Der stationare Reaktor (Punktmodell) 

Im Rahmen des sogenannten Punktmodells beschreibt man einen Kemreaktor als homogenes Objekt ohne in- 
nere Substruktur - quasi als Punkt. Im Rahmen diese Modells konnen ausgehend von einigen detaillierten 
Vorbetrachtungen gewisse weitestgehend gerechtfertigte Annahmen getroffen werden, welche in der Folge 
dargelegt werden soUen. 

In diesem Abschnitt betrachten wir vorerst den stationaren Reaktor, in welchem gewisse physikalischen Grossen 
zeitlich konstant sind - zumindest auf makroskopischen Zeitskalen der Grossenordnung Stunden oder gar 

langer. Als erste Annahme setzen wir an, dass die thermische Neutronenproduktionrate P im Reaktor durch 
Kettenreaktion in einem Reaktor proportional zur Anzahl bereits vorhandener thermischer Neutronen ist: 

P = V-n. (1) 
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Dies ist einleuchtend, da eine Veidoppelung der Neutronenzahl im Reaktor auch eine Verdoppelung der durch 
die Neutronen induzierten Kernspaltungen nach sich Ziehen wiirde, welche wiederum die Quelle fiir die Freiset- 
zung weiterer Neutronen sind. Die Produktionswahrscheinlichkeit V hat die physikalische Dimension einer 
inversen Zeit, entsprechend definiert man die sogenannte Generationszeit (mean generation time) 

A = ^. (2) 
In vielen kommerziellen Reaktoren gilt in guter Naherung A ~ 10~^s . . . 10~^s. 



AnschauUch gesprochen entspricht die Generationszeit derjenigen Zeit, in der eine gewisse Anzahl bereits 
vorhandener thermischer Neutronen in einem Reaktor mit konstant gehaltener thermischer Neutronenzahl durch 
Kernspaltung dieselbe Zahl neuer thermischer Neutronen erzeugt. Natiirlich werden bei einer (typischerweise 
durch thermische Neutronen ausgelosten) Kernspaltung zuerst schnelle Neutronen mit einer mittleren Energie 
von etwa E = 2MeV = 3.2 • 10~^^J freigesetzt. Diese werden aber innerhalb einer Zeitskala von einer 
millionstel Sekunde in den thermischen Bereich abgebremst. Thermische Neutronen werden oft durch eine 
etwas willkiirliche cwtoj^-Bedingung an ihre kinetische Energie der Art Ef^in < 5 • kBTu definiert. Dabei ist 
Tm die Moderatortemperatur und Ub = 1.38 • 10^^'^ JK^^ die Boltzmann-Konstante. 

Von den im Reaktor vorhandenen thermischen Neutronen gehen auch einige verlustig, ohne jemals eine Kern- 
spaltung auszulosen. Einerseits werden Neutronen im Reaktormaterial absorbiert, andererseits werden auch 
thermische Neutronen aus dem Reaktor hinausgestreut (Leckage). Fiir die Verlust- oder Destruktionsrate D 
setzt man mit Hilfe der Destruktionswahrscheinlichkeit V ganz allgemein an: 

D = V-n. (3) 

r = 1 /T> ist die mittlere Lebensdauer eines thermischen Neutrons im Reaktor, gemessen ab dem Zeitpunkt 
seiner Enstehung als promptes oder verz(3gert emittiertes Neutron. Stellen wir uns als Gedankenexperiment vor, 
dass in einem Reaktor plotzlich keine neuen Neutronen mehr erzeugt wiirden, so wiirde die Anzahl verbleiben- 
der Neutronen gemass der Zerfallsgleichung 

^ = -V-n n{t) = no ■ e"^'* = no • e"*/^ (4) 

abnehmen. Nach der Zeit r ware die Anzahl der Neutronen um einen Faktor e ~ 2.718... reduziert, und die 
mittlere Lebensdauer der Neutronen ist, wie aus der Theorie der radioaktiven Zerfalle bekannt, ebenso durch r 
gegeben. 



SchUesslich werden in einem Reaktor durch Zerfalle radioaktiver Substanzen fortlaufend Neutronen freige- 
setzt. Diese Substanzen entstehen beim Betrieb des Reaktors; zugleich besteht natiirlich die Moglichkeit, in 
einen Reaktor eine kiinstliche Neutronenquelle gezielt einzubringen. Die totale thermische Neutronenproduk- 
tionsrate aller vorhandenen Quellen sei im Folgenden mit Qth bezeichnet. 

In einem stationaren Kernreaktor ist die die Anderungsrate der Neutronenzahl n{t) = ^n{t) also gegeben 
durch 

h = P - D + Qth = {V -V) ■ n + Qth = . (5) 

Welter soil vorerst auch V < V gelten. Der Reaktor ist dann unterkritisch, die Kettenreaktion defizitar, da 
durch diese mehr Neutronen vernichtet werden als erzeugt. Dann kann Gleichung Q in der Form 

V -V 

■n + K-Qth = ^ (6) 

geschrieben werden. Fiihrt man die praktische Abkiirzung 

V -V 

P=^ (7) 
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ein, so gilt 

p-n + K-Qth = Q. (8) 

Tatsachlich lasst sich die Generationszeit auch naherungsweise als die typische Zeit interpretieren, die ein Neu- 
tron fiir das Durchlaufen des Fermizyklus von der Entstehung bei einer Kernspaltung bis zum Auslosen einer 
erneuten Kernspaltung benotigt. Gleichung ([8]l erhalt damit eine einfache Interpretation: Wahrend einer Gen- 
erationszeit A werden in einem Kernreaktor A • Qth Neutronen durch die Quelle(n) erzeugt, wahrend von den 
vorhandenen n Neutronen p ■ n Neutronen durch den defizitaren Fermizyklus vernichtet werden. Im stationaren 
Zustand halten sich diese beiden Prozesse die Waage, und die Neutronenquellstarke des Reaktors gleicht den 
Neutronenverlust gerade aus. 



Die Grosse p wird Reaktivitdt genannt. Offensichtlich ist p negativ, wenn V KV.lm stationaren Fall folgt aus 
Gleichung ^ 

n = ---A-Qth. (9) 

P 

Ein Reaktor mit p < heisst unterkritisch, und aus Gleichung (|9]l wird ersichtlich, dass die stationare Neutro- 
nenzahl in einem Reaktor stark ansteigt, wenn die negative Reaktivitat eines Reaktors in die Nahe des Wertes 
p = geriickt wird. An diesem Anstieg von n ~ — 1/p lasst sich messtechnisch erkennen, wann der kritische 
Zustand p = (oder P = P) eines Reaktor unmittelbar bevorsteht. 

Technisch kann die Reaktivitat beispielweise durch das Ausfahren von Steuerstaben erreicht werden. Dadurch 
werden im Reaktor weniger Neutronen absorbiert, die Destruktionsrate so gesenkt und p erhoht. Die resul- 
tierende Erhohung der Neutronenzahl im Reaktor fiihrt zur Erhohung der lokalen thermischen Neutronen- 
flussdichten, welche durch (exteme) Neutronendetektoren detektiert werden konnen. Die Neutronenfluss- 
dichte ^th = n'^f^ • v ist definiert als das Produkt der thermischen Neutronendichte mit der (mittleren) 
Geschwindigkeit v der thermischen Neutronen. 



Der Multiplikationsfaktor 

Alternativ zur Reaktivitat p ist auch die alternative formale Verwendung des Multiplikationsfaktors 

k = V/V = t/A (10) 
von praktischem Nutzen. Gemass der Definition besteht also die Beziehung 

^-P k-1 1 

P= ^ =—j—, k=- . (11) 

V k ^ — P 

Damit kann Gleichung ([8]) auch in der Form 

V 1 

- l)n = -k ■ A ■ Qth = ■ :p ■ Qth = -T ■ Qth (12) 

geschrieben werden. Innerhalb der mittleren Neutronenlebensdauer r nahme also die Neutronenzahl in einem 
unterkritischen Reaktor (mit A; < 1) von n um den Multiplikationsfaktor k auf k ■ n ab, wiirde die Quelle den 
Neutronenverlust {k — l)n nicht mit r • Qth Neutronen kompensieren. Im iiberkritischen Reaktor mit A; > 1 
wachst die Neutronenzahl auch ohne Quelle "multiplikativ" oder exponentiell an. Dieser dynamische Vorgang 
muss zum besseren Verstandnis aber einer etwas detaillierteren Betrachtung unterzogen werden. 



Die kinetischen Grundgleichungen 

Wollten wir das instationare, zeitliche Verhalten der Neutronenzahl in einem Reaktor untersuchen, so ware ein 
voreiliger Ansatz gegeben durch 

n = Y • n + QSth oder A ■ h = p ■ n + A ■ QSth ■ (13) 
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Diese Gleichung driickt folgenden Umstand aus: Wahrend einer Generationszeit werden im Kernreaktor A • h 
thermische Neutronen gebildet (re kann in dieser kurzen Zeitspanne als konstant betrachtet werden), davon wer- 
den A • Qth Neutronen durch die Quelle und p ■ re Neutronen durch die Kettenreaktion erzeugt. 



Nun werden bei einer Kemspaltung tatsachlich nicht alle Neutronen sofort, also prompt freigesetzt. Es kommt 
vor, dass bei einer Kernspaltung Bruchstiicke des urspriinglichen Spaltkerns entstehen, die erst mit einer 
zeitlichen Verzogerung in der Grossenordnung von einigen Sekunden bis Minuten schliesslich doch noch ein 
Neutron emittieren. Bereits 1939 wurde dieser Umstand bemerkt[l |. Sole he Kerne werden Vorlaufer {delayed 
neutron precursors) genannt. Bis 2001 waren 382 verschiedene Vorlaufer bekanntQ. 

In einem vereinfachten, aber zweckmassigen Modell kann angenommen werden, dass alle Vorlaufer dieselbe 
Lebensdauer besitzen. Fiir viele Betrachtungen ist eine durchschnittliche Lebensdauer von / = 13s ein prak- 
tischer Wert. Dann ist es so, dass in einem Reaktor fortlaufend Vorlaufer erzeugt werden, die wiederum Neu- 
tronen freisetzen. Experimentell findet man, dass bei einer Kernspaltung in reinem ^35^ gt^^a /3 = /3c/-235 = 
0.64% der Neutronen verzogert freigesetzt werden. 

Die resultierenden Bilanzgleichungen fiir die Neutronenzahl re und Vorlauferzahl C lauten somit 

n = + \C + Qth (14) 

und 

C=^n-\C . (15) 
A 

Diese Gleichungen werden auch (vereinfachte) kinetische Grundgleichungen der Reaktorkinetik oder reaktor- 
dynamische Punkgleichungen genannt. 



Der Term — ^ • re in der ersten Gleichung (|14|) beriicksichtigt den Umstand, dass in jedem einmal durchlaufe- 



nen neutronenzahlerhaltenden Fermizyklus der Dauer A jeweils /3 • re Neutronen fiir die Vorlauferproduktion 
abgezweigt werden. Zugleich zerfallen C momentan vorhandene Vorlaufer mit einer Rate von A • C und setzen 
dabei ein Neutron frei. Die Zerfallsrate beim radioaktiven Zerfallsgesetz ist uber die Beziehung X = l/l mit 
der mittleren Lebensdauer der zerfallenden Kerne verknupft. 



Die zweite Gleichung ( 15 I beschreibt die Produktion der Vorlaufer durch Neutronen (/? • re Vorlaufer pro Fer- 



mizyklus der Dauer A) und beriicksichtigt auch deren Zerfallsrate durch den Term — AC. 

Tatsachlich werden die prompten Neutronen mit einer durchschnittlichen Energie von etwa 2 MeV (siehe 
Anhang), die verzogerten Neutronen mit durchschnittlich etwa 0.5 MeV freigesetzt. Daher ist die Thermal- 
isierungswahrscheinlichkeit, also die Wahrscheinlichkeit, dass ein hochenergetisches Neutron in den thermis- 
chen Bereich abgebremst werden kann, fiir prompte und verzogerte Neutronen nicht wirklich gleich. Die 
Grossen C und /3 in den kinetischen Grundgleichungen entsprechen effektiven, thermisch gewichteten Grossen, 
die aber von der tatsachlichen Vorlauferzahl und dem naiv definierten /3-Wert in realistischen Fallen nicht stark 
abweichen. Eine ahnliche Aussage gilt fiir die Neutronenzahl re in den reaktorkinetischen Grundgleichun- 
gen (Punktgleichungen). Ein Neutron, welches am Rande eines Reaktors entsteht, lost auf Grund seiner ho- 
hen Leckagerate mit geringerer Wahrscheinlichkeit eine weitere Kernspaltung aus und tragi so weniger zur 
Neutronenbilanz bei als ein im Inneren des Reaktors entstandenes Neutron. Diesem Umstand kann mit Hilfe 
einer Gewichtungsfunktion oder Einflussfunktion ("Importance function") Rechnung getragen werden, aus 
welcher sich eine fiir die Punktgleichungen relevante gewichtete Neutronenzahl der Form 



re 



<S>\r) n°j^{r} dV (16) 



Reaktorkern 



errechnen lasst. Trotz diesem Umstand ist es fiir das Verstandnis der Reaktorkinetik hinreichend, wenn re als 
eigentliche Gesamtzahl der thermischen Neutronen aufgefasst wird. In guter Naherung kann <I>^ zur im Reaktor 
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herrschen Neutronenflussdichte ^th proportional gesetzt werden. 



Weiter muss bemerkt werden, dass die bei der Kemspaltung entstehenden Vorlauferkeme bei etwas praziseren 
Betrachtungen typischerweise in sechs Gruppen von Vorlaufern mit jeweils vergleichbaren Lebensdauern zusam- 
mengefasst werden. Es wird in diesem Falle die Anzahl der in Gruppe i zusammengefassten Vorlaufer, welche 
ein verzogertes Neutron freisetzen werden, durch eine effektive Anzahl Ci und eine mittlere Gruppenzerfalls- 
konstante Aj angegeben. Eine eigene "Gruppe" bildet der langlebigste Vorlaufer ®^Br mit einer Lebensdauer 
von 1/Xq = 80.387 Sekunden oder einer Halbwertszeit von 55.72s (siehe Anhang). Ein ^^Br-Kem wandelt 
sich aber mit 97%-iger Wahrscheinlichkeit durch /3^-Zerfall in einen ^''Kr-Kem um, bei lediglich 2.6% der 
Zerfalle voUzieht sich ein /3~n-Zerfall, bei dem bei einem /3-Zerfall zugleich noch ein Neutron emittiert wird - 
bei diesem Prozess entsteht stabiles ®^Kr. Die entsprechende Vorlauferzahl Cq zahlt nur diesen Anteil der fiir 
die Kettenreaktion relevanten ^^Br- Kerne. 



Im stationaren Reaktor gilt C = 0, wegen Gleichung ( 15 1 folgt 



C = = ^n-XC, esistalso C=-^n~10^n. (17) 
A AA 

Die Vorlauferzahl iibertrifft die Neutronenzahl im stationaren Reaktor bei Weitem und die Vorlaufer stellen ein 
grosses, trages Reservoir an latenten Neutronen dar. 



Offensichtlich ist im stationaren Zustand auch AC = £n und somit folgt weiter mit Gleichung ( 14i im sta- 
tionaren Zustand 

= + XC + Qth = + + Qtt = |n + Qt^ , (18) 

wie wir es in den einfiihrenden Betrachtungen dargelegt wurde. 

Zeitverhalten der Neutronenzahl im Reaktor mit konstanter Reaktivitat p 
Losung der vereinfachten kinetischen Grundgleichungen: Vorbetrachtungen 

Die vereinfachten kinetischen Grundgleichungen mit einer Vorlaufergruppe lauten filr die (gewichtete) Neutro- 
nenzahl n und die thermisch gewichtete Vorlauferzahl C bei Vemachlassigung der Quelle 

h = f^^n +XC, (19) 
A 

C = - AC, (20) 
A 

wobei A ~ 10~^ s die Generationszeit, (3 ~ 0.64% = 0.0064 den (effektiven) Anteil verzogerter (thermischer) 
Neutronen in einem mit reinem ^35^ betriebenen Reaktor und A = l/l ~ 0.0773s^^ die durchschnittliche 
Zerfallskonstante der Vorlaufer darstellt. Die Reaktivitat p kann durch externe Einfliisse kurzfristig variiert 
werden. Fiir grobe Abschatzungen kann auch A ~ 10~^ s, A = 0.1 s^^ oder ^ = 10 s und (3 = 0.01 = 10^^ 
angenommen werden. 

Bei schwach angereichertem Uran spielt die Schnellspaltung des 238^ gjj^g RoUe. Da bei der Spaltung von 
^^^U-Kemen gar etwa 1.5% der Neutronen verzogert freigesetzt werden, kann das effektive (3 in Versuchsreak- 
toren wie zum Beispiel dem CROCUS der Eidgenossischen Technischen Hochschule in Lausanne, welcher 
nur schwach angereichertes Uran enthalt, einen Wert von 0.78% erreichen. Dabei spielt auch der Umstand 
eine Rolle, dass prompte Neutronen auf Grund ihrer hohen Energie (Ekm — 2MeV) mit einer etwas grosseren 
Wahrscheinlichkeit aus einem kleinen Versuchsreaktor herausgestreut werden als die etwas energiearmeren 
verzogerten Neutronen (Efcj„ ~ O.SMeV). Da in kommerziellen Reaktoren im Verlaufe ihres Betriebs Pluto- 
nium mit /3p„_239 = 0.21% erbriitet wird, nimmt der Anteil verzogerter Neutronen im Verlaufe des Leistungs- 
betriebs wiederum leicht ab. 
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Die obige Form der kinetischen Grundgleichungen vernachlassigt ja den Einfluss einer Neutronenquelle, die 
unabhangig von der Kernspaltung mit zeitlich etwa konstanter Rate Qth thermische Neutronen im Reaktor 
erzeugt. Diese Naherung ist im iiberkritischen Reaktor (p > 0) meistens sinnvoll. Fiir p < muss Gleichung 



( [T9| ) ersetzt werden durch 

n = + AC + Qth , (21) 

mit entsprechenden Konsequenzen, die spater diskutiert werden. 

Im stationaren Zustand gilt nun n = C = 0, per definitionem sind in einem stationaren Zustand physikalisch 
messbare Grossen zeitlich unveranderlich. 

Gemass Gleichung ( [l7| gilt im stationaren Zustand C = ^n. Befindet sich ein Reaktor also im stationaren 
Zustand und wird in diesem gestort, so gilt anfanglich 

^n = XC. (22) 

Ist weiter die Reaktivitat p von ahnlicher Grossenordnung wie (3, so ist auch von ahnlicher Grossenordnung 

wie ^j^n. 

Daher andert sich die Neutronenzahl zeitlich um ahnliche Betrage wie die der Vorlaufer (solange der Reaktor 
prompt unterkritisch bleibt: p < P beziehungsweise pp = p — P < 0) doch diese sind etwa tausendfach 
zahlreicher. Entsprechend reagiert die Vorlauferzahl etwa tausendfach langsamer als die Neutronenzahl auf 
Storungen des stationaren Zustandes. Diese Aussage gilt allerdings nicht mehr im prompt iiberkritischen Reak- 
tor mit p > (3,m dem die Neutronenzahl schon kurze Zeit nach dem Erreichen der gefahrlich hohen Reaktivitat 
die Vorlauferzahl iibertrifft. 

Die entsprechenden Zeitskalen im prompt unterkritischen Reaktor sind fiir den Bereich < p < 0.9/3 gegeben 
durch die prompte Zeitkonstante copr 



te,pr = = ^ ^7T7T7 =0-01---0-ls, (23) 



10-S 
P- p~ 0.01... 0.001 

dabei nennt man te^pr auch die prompte Zeitverzogerung, die sogenannte verzogerte Zeitkonstante ist gegeben 
durch ujy = = A ~ 0.1 s~^. tf,^v heisst verzogerte Zeitverzogerung, auf dieser Zeitskala reagieren die 

Vorlaufer in der Eingruppennaherung auf Anderungen der Neutronenbilanz. 

Bevor die exakten Losungen der Eingruppengleichungen untersucht werden, soil an dieser Stelle kurz die Situa- 
tion diskutiert werden, bei der die Reaktivitat po in einem Reaktor mit no Neutronen abrupt einen neuen Wert 
p annimmt. Vereinfachend soli dabei gelten, dass innerhalb des betrachteten Zeitrahmens die Vorlauferzahl als 
eine Konstante betrachtet werden darf. Aus 

n = ^^ra + AC = + AC (24) 
mit der prompten Reaktivitat Pp < folgt so 

n(t)=fno- — -cV^^ + ^-C. (25) 

\ pp J Pp 



Dieser zeitliche Verlauf der Neutronenzahl lost offensichtlich Gleichung ( 24 1 mit der Anfangswertbedingung 
n{t = 0) = no- Die Losung zerfallt mit der charakteristischen Zeitskala gegen den (quasi)stationaren 
Gleichgewichtswert n{t » te,pr) = {^^/ Pp) ■ C. Letztlich reagiert aber auf einer grosseren Zeitskala auch die 
Vorlauferzahl auf die veranderten Bedingungen. 
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Losung (mit spezieller Betrachtung des verzogert kritischen Zustands < p < /3) 



Um die vereinfachten kinetischen Grundgleichungen exakt zu losen, leiten wir Gleichung ( 19 1 zuerst noch 
einmal ab und erhalten flugs 



n 



A 



Aus der Summe von Gleichung (19) und (20 1 folgt aber sofort 



C + h = —n Oder AC 
A 



\p 
A 



n — Xh , 



ersetzen wir also AC in Gleichung (26 1 durch den Ausdruck aus Gleichung (27 ), so ergibt sich 

.., /3 -/>.,> . A/3 

n H : — n + An —n = U . 

A A 



(26) 



(27) 



(28) 



Diese Gleichung lasst sich durch einen Exponentialansatz einfach losen - wir setzen n{t) = uq ■ e'^ * und 
erhalten wegen ^e"^ * — - • 



u! ■ e 



n = oj^ ■ n = oj^ ■ riQ ■ e^'^ und n 



u ■ n = uj ■ uq ■ e 



Dies eingesetzt in Gleichung (28 1 ergibt 



u? -riQ- e'^'* + bj- 



P-P 
A 



+ A -no • e' 



(29) 



(30) 



Natiirlich konnen die Faktoren no und e^ * aus dieser Gleichung weggekiirzt werden. Somit erhalten wir die 
sogenannte Sakulargleichung fiir die reziproke stabile Reaktorperiode (Zeitkonstante) 



A 



+ X\uj 



Xp 
A 



0. 



(31) 



Aus der quadratischen Gleichung (31 1 ergeben sich zwei mogliche Losungen fiir uj, namlich 

l( P-P . \ . 1 



A 



A ±: 



p-l 

A 



X +4 



A 



(32) 



Folgende Eigenschaft dieses Audruck ist bemerkenswert, sofern die Reaktivitat des Reaktors positiv ist - diesen 
Fall wollen wir in der Folge zuerst diskutieren. Der Wurzelterm ist namlich betragsmassig grosser als der linke 
Klammerterm in Gleichung ( 32 1, da unter der Wurzel zum positiven Quadrat des Klammerterms auch noch eine 
positive Grosse 4Ap/A hinzu addiert wird. Daher ist ioi (definiert mit dem Pluszeichen in der Losungsformel) 
immer positiv, und uj2 ( "-") immer negativ. Die beiden Losungskonstanten entsprechen also einer exponen- 
tiell ansteigenden und einer exponentiell abfallenden Komponente in der Neutronenzahl, die sich im Reaktor 
uberlagern. Die exponentiell gedampfte Losung kann aber bei stabiler Reaktorperiode vernachlassigt werden, 
wenn also im iiberkritischen Reaktor die Neutronenzahl exponenziell zu wachsen beginnt. Offensichtlich ist 
U2 mit der prompten Zeitverzogerung eng verkniipft; sie wird auch als transiente Reaktorperiode bezeichnet. 
Die allgemeinste Losung der Gleichung ([28]) lautet also 



n{t) = Uex ■ e"'!-* + n • e'^^- 



(33) 



wobei der e'^^ ^-Anteil exponentiell abgedampft ist. 

Weiter ist im Bereich < p < 0.9 • /3 (oder < /9„ < 90^ 

bedeutend grosser als A. Man kann daher vereinfachen: 



0,9$) der betragsmassige Wert von 



1 



\ 



p-p 

A 



A +4 



Xp I p-p 



A 



A 



+ A 



'1 + 



4Xp A2 



(34) 
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Da (^pj^p^i viel kleiner als eins ist und t/I + x ~ 1 + ^2; fiir kleine x, folgt schliesslich 
(wieder mit A < (/? - p) /A) 



2 V 2 \ A 
Somit ist uj = (jJi naherungsweise durch die Formel 



gegeben. 



(35) 



(36) 



Die stabile Reaktorperiode ist also im Rahmen der vereinfachten kinetischen Grundgleichungen naherungsweise 
gegeben durch 



te 



Xp ' 



die exakte Formel ist aber gemass Gleichung ( 32 1 



\ 



Die stabile Reaktorperiode ist gleichwertig zur reziproken stabilen Reaktorperiode uj 

1 



te 



UJ 



(37) 



(38) 



(39) 



die Verdopp(e)lungszeit, d.h. die Zeit, nach der sich im eingeschwungenen Zustand die Neutronenzahl jeweils 
verdoppelt, ist etwas kleiner 

(40) 



t2 = ln(2) • te = —-^ ~ 0.69 • te . 



UJ 



Figur 2 stellt die Beziehung aus Gleichung ([38]) graphisch dar. Im angezeigten Bereich lassen sich die Unter- 



schiede zwischen Gleichung (37 1 und (38 1 von Auge kaum ausmachen. 



Inhourbeziehung (Eingruppengleichung) 




Figur 2: Verdoppelungszeit versus Reaktivitat in der Eingruppennaherung. 
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Die Inhourbeziehung 

Betrachten wir die Inhourkurve (Figur 3), welche den Zusammenhang zwischen Reaktivitat und Verdop- 
pelungszeit eines typischen Reaktors mit /3 = 0.64% darstellt, so fallt primar eine gewisse Diskrepanz im 
Bereich kleinerer Verdoppelungszeiten gegeniiber Figur 2 auf. Der Grund dafiir liegt in der simplen Tatsache 
begriindet, dass eine prazisere Beiechnung der Verdoppelungszeit die Beriicksichtigung aller 6 Vorlaufergruppen 
mit ihren individuellen Zerfallkonstanten voraussetzt. Die voUen kinetischen Grundgleichungen 



n = — . — n 
A 



+ Y,X^Ci (41) 



C^ = ^n- Xid (42) 

lassen sich tatsachlich durch einen Exponentialansatz losen. Es stellt sich heraus, dass sich die Neutronenzahl 
als Summe von sieben verschiedenen Exponentialfunktionen 

7 

n{t) = riex ■ e"^^'* + ^nj ■ e'^^'* (43) 
i=2 

schreiben lasst. Eine analoge Aussage gilt fur die Vorlauferzahlen. Dabei folgt aus einer detaillierteren mathe- 
matischen Analyse, dass im Falle p > 

w = a;i > und < fur j = 2, . . . 7 . (44) 

Fiir /9 < sind alle ujj negativ, im unterkritischen Reaktor zerfallt jede Storung der Neutronenzahl, und so strebt 
diese jeweils gegen den durch die Neutronenquellen vorgegebenen Gleichgewichtswert gemass Gleichung (|9]). 
Fiir p > existiert aber eine stabile Reaktorperiode welche in der Inhourkurve (Figur 3) abgelesen werden 
kann. Die restlichen Exponentialfunktionen fallen zeitlich stark ab und beschreiben den Einschwingvorgang, 
der einem stabilen exponentiellen Anstieg der Neutronenzahl zeitlich vorangeht. Nach dem Einschwingvorgang 
wird die Entwicklung der Neutronenzahl faktisch nur noch durch den Term 

n{t) = Hex ■ e"^-* (45) 

7 

beschrieben. Offensichtlich ist n(0) = n^x + 1^ fij die Startzahl der Neutronen, wobei natiirlich bei unseren 

J=2 

Betrachtungen immer p =konst. fiir t > gilt. In einem Leistungsreaktor ist es typischerweise so, dass eine 
Reaktivitatserhohung aus dem kritischen Zustand p = zu p > eine Erhitzung des Brennstoffs (und des 
Kiihlmittels oder Moderators) bewirkt, welche die Reaktivitat des Reaktors wieder reduziert. Durch diese sta- 
bilisierende Ruckkopplung geht der Reaktor bei erhohter Leistung letztlich wieder in den kritischen Zustand 
iiber. 



Die Zeitkonstante Lj02 im Eingruppenmodell kann wegen A ^ {13 — p)/ Kin Gleichung (32) durch 



U2 ^ (46) 

angenahert werden. Es handelt sich bei ui2 (bis auf ein Vorzeichen) namlich um die prompte Zeitkonstante 
ujpr- Sie reprasentiert im Wesentlichen die zusatzlichen 6 Zeitkonstanten a;2, . . . , des voUen Modells mit 6 
Vorlaufergruppen. Diese Vereinfachung wird in gewissen Fallen einer realistischen Beschreibung nicht mehr 
gerecht. 

Man kann zeigen, dass die Reaktivitat des Reaktors mit jeder Zeitkonstanten Wj fur j = 1, 7 die exakte 
mathematische Beziehung 

^A + ^^'j (47) 
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erfiillt. Kann in dieser Gleichung A vernachlassigt werden, so erhalten wir fiir die Eingruppengleichungen 



p = iO 



0J + \ 



A, 



(48) 



das heisst, wir gelangen wieder zu Gleichung (37i, wobei die negative reziproke Periode in der gemachten 
Naherung unter den Tisch fallt. 



Die stabile Reaktorperiode ist also naherungsweise durch 

p \ p 

gegeben, und die Verdoppelungszeit folgt aus der Halbwertszeit = • ^ — 8, 97s der Vorlaufer 



(49) 



t2 



■ t 



P 



1/2 



^ P = P- 



t 



1/2 



t2 + t 



1/2 



(50) 



Misst man also die stabile Neutronenzunahme in einem Reaktor durch Messung eines Neutronenflusses $ cx n, 
so lasst sich also mit Hilfe von Formel ( 50 1 sofort die Reaktivitat naherungsweise berechnen. Genauere Aus- 
sagen gewinnt man aus der Inhourkurve, welche alle Vorlaufer berilcksichtigt. 

Ein Beispiel: Fiir eine Reaktivitat von 80pcm erhalten wir 

0.0064 - 0.0008 



t2 



0.0008 

in passabler Ubereinstimmung mit t2 = 50s in der Inhourkurve. 



9s = 63s , 



(51) 




20 30 40 50 100 200 t2[s] 

Figur 3: Die Inhourkurve. 

Der Ausdruck Inhour geht nicht aus einen Forscher gleichen Namens zuriick, sondern beruht historisch und 
vereinfacht gesagt auf der Tatsache, dass bei friihen Reaktorexperimenten (— )^ CP — 1) der Neutronenzuwachs 
pro Stunde, also in Beziehung zu inverse hours, betrachtet wurde. 

Da sich die Neutronenzahl im Bereich < p < /? in zeitlich in kontroUierbarer Weise vergrossert, nennt 
man den entsprechenden Reaktivitatsbereich verzogert kritisch. Es ist also der Existenz der Vorlaufer zu ver- 
danken, dass sich die Leistung eines Reaktors zeitlich im Minutenbereich regeln lasst, sofem in einem sicheren 
Reaktivitatsbereich unterhalb von ~ 200pcm operiert wird. 
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AUgemeine Diskussion der Losungen der kinetischen Grundgleichungen 
Die Reaktivitat als Ordnungsparameter 

In der Folge soUen nach der eingehenden Diskussion des iiberkritischen Reaktors nun auch noch die weiteren 
durch die Reaktivitat charakterisierbaren Reaktorzustande diskutiert werden. Ein Reaktor mit negativer Reak- 
tivitat p < heisst also unterkritisch, bei p = ist er kritisch; gilt /> > 0, so ist der Reaktor uberkritisch. 
Bei iiberkritischen Reaktoren unterscheidet man den kontrollierbaren Fall des Reaktors im Anfahrbereich mit 
< p < /3, welcher auch als verzogert uberkritisch bezeichnet wird. Der prompt kritische (p = (3) oder prompt 
iiberkritische Zustand {p > /3) muss auf jeden Fall vermieden werden. 

Quellbereich (p < 0) 

Im unterkritischen Bereich sind die beiden inversen Reaktorperioden 

^i,2 = 2(^^-Al±^, l^-A) +4:f . (52) 

negativ. Ohne Neutronenquelle wiirde also die Neutronenzahl in einem Reaktor gegen Null zerfallen. Beriick- 
sichtigt man eine potenziell vorhandene Neutronenquelle, so kann aus der allgemeinsten Losung der verein- 
fachten kinetischen Grundgleichungen ohne Quelle 

n{t) = m ■ + na • e"^-* (53) 

sofort eine Losung des Gleichungssystems 

n = +\C + Qth , (54) 





C=^n-XC (55) 
^^LL^ + „,.e-i-* + „2.e--* , (56) 



angegeben werden: 



C(t) = -f:-Qi,+ ( ^ + •ni.e--*+ l^^ + ^j -na-e--* . (57) 

Fiir grosse Zeiten zerfallt diese Losung zu den stationaren Werten 

_ K-Qth ^ _ P n - ^ 
p Xp AA 

welche sich sofort aus den stationaren, zeitunabhangigen Bilanzgleichungen 

n = = ^^^^ +\C + Qth, C = Q=^n-\C (59) 

ergeben. Offensichtlich ist die stationare Neutronenzahl und die Vorlauferzahl im Reaktor proportional zu 
— ^ . Entsprechend dauert es bei Reaktivitatsanderungen im unterkritischen Reaktor umso langer, bis sich der 
Reaktor wieder auf einen stationaren Zustand eingeschwungen hat, je naher man sich an der kritischen Schwelle 
p = befindet. Wird die negative Reaktivitat zu betragsmassig immer kleineren Werten reduziert, so stellt sich 
nach einiger Zeit eine immer grosser werdende stationare Neutronen- und Vorlauferzahl ein. Der Aufbau 
dieser zunehmend grosser werdenden Zahlen dauert zunehmend langer, da die eigentliche Kettenreaktion im 
unterkritischen Reaktor immer defizitar ist - das gekoppelte System der Neutronen und Vorlaufer wird letztlich 
aus der Reaktorquelle erzeugt, die mit konstanter Rate Neutronen nachliefert. Bildet man namlich die Summe 



der Gleichungen ([54]) und (55 1, so ergibt sich 



ri + C = ^n + Qth. (60) 
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Der Term -^n ist aber negativ, da p negativ ist. Die Gesamtzahl ii + C der Neutronen und Vorlaufer zusammen, 



die f iir schwindendes p immer grossere Werte annimmt, kann nur dank des Quellterms Qth in Gleichung ( 60 1 
wachsen. 



Bei einer Reaktivitatsanderung im Bereich betragsmassig kleiner Reaktivitaten ergibt sich die Einschwingzeit 
tes eines Reaktors im Wesendichen aus dem am langsamsten zerfallenden Exponentialterm in Gleichung ( 56 1 



(siehe auch Gleichung (36 1) mit reziproker Periode 



> 





A +4 



Xp 
A 



Ap 

P-p 



Xp 



(0 < -p « /?) 



zu 



1 

UJl 



A. 

A|H 



(61) 



(62) 



I . Der tatsachUche Ein- 



Fiir /? = 0.0064 = 640pcm und 1/A = 12.9s ergibt sich daraus t^s — |p| inpcm 
schwingvorgang dauert natiirlich (5-lOmal) langer als t^s, da die Eingruppengleichung den Einfluss der lan- 
glebigen Vorlaufer unterschatzt und da tes im Wesentlichen die Zeit fiir eine Reduktion der Abweichung des 
Neutronen- Vorlaufer-Systems vom Gleichgewichtszustand auf einen e-tel beschreibt. 

Wird die Reaktivitatsanderung schnell (prompt) durchgefiihrt, so vollzieht die Neutronenzahl zuerst einen soge- 
nannten prompten Sprung, bei welchem ein quasistationares Gleichgewicht zwischen Neutronen und Vorlaufer 
herbeigefiihrt wird. Dieser Vorgang wird mit der betragsmassig grosseren Zeitkonstante 0J2 zeitlich beschrieben. 
Bei kleinen Reaktivitasanderung ist der prompte Sprung weniger ausgepragt. Quantitative Betrachtungen zum 
prompten Sprung finden sich im Abschnitt Schnelle Reaktivitdtsdnderungen. 



Kritischer Reaktor: p = 

Im kritischen Reaktor vereinfachen sich die kinetischen Grundgleichungen mit p = zu 



n 



A 



n +\C + Q 



th, C = ^n- XC 



und bildet man die Summe dieser beiden Gleichungen, so folgt 

C + h = Qth, 

also 

C + n = No + Qth-t. 



(63) 



(64) 



(65) 



Da die Gesamtzahl von Vorlaufern und Neutronen offensichtlich linear mit der Zeit ansteigt, da die Kettenreak- 
tion eigentlich selbsterhaltend ist, die Neutronenquelle aber immer mehr Neutronen in den Reaktor einspeist, 
liegt es nahe, sowohl fiir die Anzahl der Vorlaufer sowie fur die Anzahl der Neutronen einen zeitlich linear 
ansteigenden Verlauf anzusetzen: 



nit) = no + Qi-t, C{t) = Co + Q2-t, Qi + Q2 = Qth 
Eingesetzt in die Neutronenbilanzgleichung ([63]) liefert der Ansatz 



Qi 



/3 (3 
-no 



Qi-t + XCo + XQ2-t + Q 



X " A 

Die in t linearen Terme miissen sich offensichtlich aufheben, also gilt 



th ■ 



Qi = XQ2 



Qi 



Q2 



XA 



Qi — Qth — Qi , 



_Q_ 
1 + 



th 



AA 



Q2 



_Q 
1 + 



th 
AA 



(66) 
(67) 

(68) 
(69) 
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In einem kritischen Reaktor mit Quelle nimmt die Neutronenzahl nach dem Einschwingen linear zu, und zwar 
in einem Zahlenverhaltnis zu den Vorlaufern, welches bereits vom stationaren, unterkritischen Fall her bekannt 
ist. Die Anfangswerte no und Co konnen geeignet gewahlt werden, sodass die kinetischen Grundgleichungen 
erfiillt werden. 

Beachte: Zur linearen Losung der kinetischen Grundgleichungen lasst sich noch eine exponentiell zerfallende 
Losung ~ e"'^*'"'* * mit inverser Periode ujf^rit = /3/A + A ~ A hinzuaddieren. Der Wert dieser Frequenz 
folgt wiederum aus Gleichung ( [32| ) mit Minuszeichen und ist auch fur den sogenannten inhomogenen Fall der 
kinetischen Grundgleichungen mit Neutronenquelle gultig. 

Der prompt kritische Reaktor (p = /3) 

Die inversen Reaktorperioden lassen sich durch 

A /A2 A/3 M 

approximieren. Die positive Reaktorperiode beschreibt einen starken Anstieg der Neutronenzahl auf der prompt 
kritischen Zeitskala tpk = — 0-2s. Da auf dieser kurzen Zeitskala die RoUe der langlebigen Vorlaufer 

durch die mittlere Zerfallskonstante A ~ O.ls^^ ubergewichtet wird, muss tpk mit einer effektiven Zerfallskon- 
stanten A ~ 0.407s^^ berechnet werden. 

Der prompt iiberkritische Reaktor (p > /3) 

Im prompt iiberkritischen Reaktor kann der Einfluss der Vorlaufer eigentlich vernachlassigt werden, da die 
prompten Neutronen allein sich ausserst schnell vervielfachen. Ist p hinreichend gross, sodass gilt 

f = ^»A, (71) 

so ist 

Bereits fiir /? = 2/3 oder pp = /3 ist l/tJi ~ 0.016s, d.h. in dieser Zeit ver-e-facht sich die Neutronenzahl im 



Reaktor. Dieses Resultat ergibt sich auch aus der Neutronenbilanzgleichung ( 19 1, wenn der Vorlauferterm AC 
vernachlassigt wird: 

n = —7- n{t) = uq ■ e~^'^ . (73) 

Natiirlich ist der prompt (uber)kritische Zustand eines Reaktors in jedem Falle zu vermeiden, da die entsprechende 
Leistungsexkursion zur Kemschmelze fiihren kann. Auf der obigen sehr kurzen Zeitskala ist es allerdings nicht 
mehr statthaft, von einer Gesamtreaktivitat (lumped reactivity) des Reaktors zu sprechen (siehe Abschnitt Dif- 
fusion). 

Schnelle Reaktivitatsanderungen 

Bei einer schnellen Reaktivitatsanderung (po = fiir t < — )• pi = p fiir t > 0) in einem kritischen Reaktor 
(zum Beispiel durch Stabeinwurf) wird die Neutronenzahl no fiir t < abrupt zu ni fiir t > verandert, 
wahrend das trage Reservoir der Vorlaufer von der Reaktivitatsanderung vorerst wenig beeinflusst wird. Innert 
kurzester Zeit stellt sich wieder ein quasistationares Gleichgewicht zwischen den Neutronen und den Vorlaufern 
ein, und es gilt in guter Naherung wie im stationaren Zustand 

— — m = ^no = AC = — ni = ^^^^i , (74) 

also ist in guter Naherung 

ni = no . (75) 

/3-p 
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Suggestiver kann man fiir den prompten Sprung mit den prompten Reaktivitaten auch schreiben {pp = p — (3) 



Pp,o ■ no = Pp,i ■ ni . (76) 
In Rahmen der Eingruppennaherung wird der zeidiche Verlauf der Neutronenzahl durch 



n{t) = no 



B _A^t p 13-p. 



p-p 



ill) 



beschrieben. Esgiltn(t = 0) = no- Da der zweite Exponentialterm sehr schnell zerfallt - die Zeitkonstante -^^^ 
(die prompte Zeitverzogerung) ist im verzogert kritischen Bereich von der Grossenordnung von Bruchteilen 
einer Sekunde - schwingt sich die Neutronenzahl sehr schnell auf 

n(f) ~ (78) 

ein. Der erste Exponentialterm beschreibt das kontroUierte Anwachsen der Neutronenzahl im verzogert kri- 
tischen Reaktor mit stabiler Reaktorperiode. Allerdings stellt sich bei realistischer Beriicksichtigung aller 
Vorlaufer die stabile Reaktorperiode erst nach einem Einschwingvorgang ein, dessen Zeitskala mit den ver- 
schiedenen Lebensdauern der Vorlaufer verkniipft ist. 

Reaktorschnellabschaltung 

Befindet sich ein Reaktor vor einer Schnellabschaltung im kritischen stationaren Zustand, so lasst sich die 



Anzahl der Vorlaufer in der i-ten Vorlaufergruppe mit Hilfe der Vorlauferbilanzgleichung (20 1 sofort angeben: 



C = ^ C, = -^n. (79) 

Durch die starke Reaktivitatserniedrigung fallt die urspriinglich vorhandene Neutronenzahl no prompt auf den 
mit den Vorlaufern im quasistationaren Gleichgewicht befindlichen Wert 

ni = no . (80) 

13-p 

Anschliessend zerfallen die Vorlaufer mit verschiedenen Zerfallskonstanten und setzen noch Neutronen frei. 
Deren Anzahl berechnet sich zu (n(0) = ni) 

In der obigen Formel wird der geringe Effekt der emeuten Vorlaufererzeugung durch die Freisetzung von 
Neutronen und nachfolgender Kernspaltung durch die noch vorhandenen Vorlaufer vernachlassigt. In guter 
Naherung gilt in Gleichung (43 1 entsprechend ujj = — Aj+i. 



Der langlebigste Vorlaufer ist ^''Br. Nachdem die kurzlebigeren Vorlaufer zerfallen sind, dominiert ^^Br die 
obige Summe und nach etwa einer Minute gilt 

n(t) ~ — • 66 • e~^« * mit 66 = 0.032, A6 = l/(80.387s) . (82) 

p p 

Dieser anfanglich exponentielle Zerfall der Neutronenzahl stabilisiert sich schliesslich nach einigen Minuten 
bei dem endlichen Wert, der durch die Quellstarke der Reaktors vorgegeben wird. 
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Die Sechsfaktorformel 



Die Sechsfaktorformel geht in vereinfachter Form als Vierfaktorformel auf Eugene Wigner0| und Enrico 
Fermi Bll zuriick und veranschaulicht die Bedeutung des Multiplikationsfaktors k aus Gleichung ([TO]). Die 



Formel beschreibt die Anzahl der Neutronen, die durch ein im Brennstoff eines Reaktors absorbiertes Neutron 
erzeugt werden und wieder den Weg zuriick in den Brennstoff finden, wo sie schliesslich erneut absorbiert 
werden. Die Formel lautet 

k = r]-e-Ps-p-Pth-f- (83) 

Dabei beschreibt r] die Neutronenergiebigkeit des Brennstoffs. Nicht jedes Neutron, welches im Brennstoff 
eines Reaktors absorbiert wird, lost namlich eine Spaltung aus. Uber ein Fiinftel der thermischen Neutronen 
wird in angereichertem Uran (3% ^^^U, 97% ^'^'^U) eingefangen, ohne eine Spaltung auszulosen. Dabei entste- 
hen die Isotope ^36^ ^j^^ ^'^^U, wobei das ^^^U schliesslich zu ^^^Pu zerfallt. Somit setzt ein im Brennstoff 
eines Druck- oder Siedewasserreaktors absorbiertes Neutron im Mittel etwas weniger als zwei Neutronen frei, 
obwohl bei einer Kernspaltung etwa i/ ~ 2.5 Neutronen entstehen. 

Die bei einer Kernspaltung entstehenden hochenergetische Neutronen konnen mit geringer Wahrscheinlichkeit 
auch ^^®U-Kerne spalten; dadurch erhoht sich die Anzahl der bereits vorhandenen Neutronen um einen Schnell- 
spaltfaktor e. Zugleich entkommen einige der schnellen Neutronen mit der schnellen Entkommwahrschein- 
lichkeit dem Reaktor und gehen verloren. Der schnelle Verbleibfaktor Pg tragt diesem Umstand Rechnung. 

Beim anschliessenden Moderationsprozess verlieren die Neutronen durch Stosse (beispielsweise an Moderator- 
kernen wie ^^C, oder ^H) an kinetischer Energie. Im Energiebereich von 7eV...200eV laufen sie dabei 
Gefahr, durch sogenannte Resonanzabsorption durch ^'^^U-Kerne eingefangen zu werden. Es kommt nur ein 
gewisser Anteil der Neutronen im thermischen Bereich an, dieser ist gegeben durch die Resonanzentkomm- 
wahrscheinlichkeit p. 

Die thermischen Neutronen diffundieren nun durch den Reaktorkem. Der thennische Verbleibfaktor Pth 
beschreibt den Anteil der Neutronen, die nicht aus dem Reaktorkem hinausdiffundieren. Die thermische 
Nutzung schliesslich beschreibt den Anteil der thermischen Neutronen, welche schliesslich wieder im Brennstoff 
und nicht im restlichen Material der Reaktorkems absorbiert werden. 

Realistische Werte fiir einen kritischen Leistungsreaktor bei konstanter VoUleistung waren in obiger Reihen- 
folge 

k = l = 1.69 • 1.04 • 0.96 • 0.787 • 0.996 • 0.756 . (84) 
Offensichtlich vergrossert sich die Neutronenzahl in einem uberkritischen Reaktor beim Durchlaufen eines 



oben beschriebenen "Fermizyklus" um den Faktor k, den wir mit k aus Gleichung ( lOl identifizieren woUen. 



Man mag sich nun fragen, wie lange ein solcher "Fermizyklus " dauert. Im Rahmen der Eingruppengleichung 
berachten wir hierzu einen schwach uberkritischen Reaktor der Reaktivitat /9 ^ /3, in welchem die Neutronen- 
zahl mit stabiler Reaktorperiode zeitlich exponentiell wachst. Dann gilt in guter Naherung fiir die reziproke 



Reaktorperiode gemass Gleichung (36 1 



- (85) 



Innerhalb der Zeit t p soil sich also die Neutronenzahl um den Faktor k erhohen. Mit einer stabilen Reaktor- 
periode gilt also wegen e^~l + xfiir|x|^l 

fc = e'^-*-^ ~ 1 + a; -tp = ~ 1 + p, (86) 

I- p 

also ist 

ti. ~ ^ ~ ^ = /3 • Z . (87) 



In Gleichung (87i taucht die mittlere Lebensdauer / ~ 13s der Vorlaufer auf. Diese dominiert die mittlere 



Dauer, die ein Fermizyklus in Anspruch nimmt. Zwar durchlauft die Mehrzahl der Neutronen den Fermizyklus 
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innert kurzer Zeit der Grossenordnung der Generationszeit A. Der geringe Anteil von weniger als einem Prozent 
der Neutronen, welche in Vorlaufern gespeichert den Fermizyklus verzogert, lebt aber ungleich langer und 
bestimmt im verzogert iiberkritischen Reaktor die mittlere Generationsdauer der Neutronen, bei der sich die 
Neutronenzahl um den Multiplikationsfaktor k erhoht. 

Diffusion 

Die Neutronenflussdichte (^th in einem Reaktorkem ist eine raumlich variierende Grosse. Typische Werte der 
Neutronenflussdichten in Leistungreaktoren bewegen sich im Bereich von einigen Raumliche 
Gradienten der Flussdichte V^^th fiihren zum lokalen Transport thermischer Neutronen, welcher durch die 
Neutronenstromdichte jth{r) im Rahmen der Diffusionsgleichung (v sei konstant) 

jth = -D ■ V^th = -D ■ grad <^th = -D-v- gradn^^ (88) 

beschrieben werden kann. D spielt hierbei die Rolle einer Diffusionskonstanten mit der Dimension einer Lange. 
Der Neutronendiffusionsstrom trachtet nach einer Dichteausgleich durch den Transport von Neutronen aus Re- 
gionen hoherer Dichte in Regionen niedrigerer Dichte. Die lokalen Neutronendichteschwankungen durch reine 
Transportphanomene (d.h. bei unrealistischer Vernachlassigung der Neutronenproduktion durch fortwahrende 
Kernspaltungen) konnen mit Hilfe der Kontinuitatsgleichung 

dfi^ —* -t -* 

-^ = -^-jth = -divjth (89) 
dt 



Oder 



beschrieben werden. Aus der Kombination von Gleichung (88 1 und (89 1 ergibt sich die Diffusionsgleichung 
(heat equation) 

-T^ = D- A^th (90) 
ot 

^-An^, D = D.v. (91) 

Die Greensche Funktion (oder Kern, Fundamentallosung) dieser Gleichung im homogenen Raum mit konstan- 
ter Diffusionskonstante lautet 

Sie beschreibt das diffusive Zerfliessen einer zur Zeit t = punktformigen Storung der Neutronendichte. Die 

I 1 d2 

Zeitskala, mit der die Gausssche Glockenfunktion in Gleichung ( 92 1 zerfliesst, ist gegeben durch — I75 > 



wenn R eine charakteristische Ausdehnung der betrachteten Konfiguration beschreibt. Fiir D ~ 1cm, R ~ 4m 
und V = 3000m/s ergibt sich to — 0.1s - damit gleichen sich Schwankungen der Neutronendichte in einem 
prompt unterkritischen Reaktor relativ schnell gegeniiber der stabilen Reaktorperiode aus, doch im prompt 
(uber)kritischen Reaktor hat das Punktmodell seine Giiltigkeitsgrenzen (auch aus anderen Griinden) erreicht. 

Anhang 

Daten der sechs Vorlaufergruppen von ^^^U 

Die Daten fiir 235^ stammen aus einer Arbeit von Keepin, Wimett und ZeiglerQ. 

Zur Illustration: Es ist tl^^f = 16.3s, tJ'J = 24.2s. Es gih Pi = P ^ 0.0064, 1 = Pik ^ 12.93s. 

i=l 1=1 

Mittlere Neutronenenergien Ei der i-ten Vorlaufergruppe sind E2 = 0.42MeV, £'3 = 0.62MeV, E^^ = 0.43MeV, 
E^ = 0.56MeV und Ea = 0.25MeV. 

Die Anteile verzogerter Neutronen fiir weitere Spaltkeme sind gegeben durch: 
Pu-2ZH = 1.57%, /3r/x-232 = 2.2%, /3[/_233 = 0.27%, PPu-2m = 0.21%. 
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Gruppe 


Vorlaufernuklide 


Relativanteil 

Pi/ P 


Zerfallskonstante Aj 
L» J 


Lebensdauer li 


Halbwertszeit ti/2,i 


1 


93Br, ... 


0.042 


3.0137 


0.332 


0.230 


2 


1401, 145Cs... 


0.116 


1.1360 


0.880 


0.610 


3 


i44cs, 1391, 90Br, ... 


0.396 


0.3014 


3.318 


2.300 


4 


1381, 89Br, ... 


0.195 


0.1114 


8.974 


6.222 


5 


i37l, 88Br, ... 


0.219 


0.03051 


32.778 


22.72 






U.UJz 


A A1 0/1 /I 

U.U1Z44 


OA ^C? 


JJ. /z 



Das Spektrum der prompten Spaltneutronen 

Nach Watt||6| lasst sich das kinetische Energiespektrum der bei der Kemspaltung prompt freigesetzten Neutro- 
nen durch den Ausdruck 

X{E) = 0.484 sinh V2Ee-^ (93) 

parametrisieren. Dabei wird die Neutronenenergie E in MeV gemessen. In ein infinitesimales Energieintervall 
werden dann je Spaltung I'iE) • x{E)dE Neutronen erzeugt, wobei die Neutronenausbeute z/ selbst von der 
Energie der Neutronen abhangt, welche die Spaltung auslosen. Neutronen, die mit hoher Energie eine Spaltung 
auslosen, setzen auch eine grossere Anzahl Spaltneutronen frei. Dieser Zusammenhang kann durch eine Uneare 
Beziehung 

u{E) = vo + a-E (94) 

dargestellt werden und wurde experimentell im Energiebereich bis 15MeV untersucht (fT'l, siehe nachfol- 
gende Tabelle). In guter Naherung lasst sich das Spaltneutronenspektrum auch durch eine Maxwell- Verteilung 
darstellen: 



:E^'\-^/^ . 



(95) 



7rV2r3/2 

mit T = 1.29MeV fiir ^^Sy. Eine approximative Maxwellverteilung der emittierten Neutronen lasst die Vorstel- 
lung eines Nukleonengases in einem hochangeregten Kern zu, der bei einer Temperatur von 1.5 • IQi^K Neu- 
tronen abdampft. Die Verteilung x^^^{E) in Gleichung (95 1 ist gemass 



x''^iE)dE=l 



(96) 







normiert. 

Der Parameter T = 1.29MeV ist mit der mittleren Energie der Spaltneutronen iiber die Beziehung 



E 



Ex''^{E)dE = 



(97) 



verkniipft. Die mittleren Neutronenenergien sind fiir die relevanten Spaltstoffe im Wesentlichen gleich: Fiir 
obiges T ist Eu-235 = 1-935 ± 0.05MeV, weiter ist Eu-233 = 1-96 ± 0.05MeV und Epu-239 = 2.00 ± 
0.05MeVL7J. 
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Kernbrennstoff 


1^0 


a 

[MeV^^] 


Energiebereich 
[MeV] 


232Th 


1.87 


0.164 


0<E<15 


233u 


2.48 


0.075 


0<E<1 




2.41 


0.136 


1<E<15 


235u 


2.43 


0.065 


0<E<1 




2.35 


0.150 


1<E<15 


238^ 


2.30 


0.160 


0<E<15 


239PU 


2.87 


0.138 


0<E<1 




2.91 


0.133 


1<E<15 



Spektrum prompter Spaltneutronen 
0.4| ^ < < 




2 3 4 

Neutronenenergie (MeV) 



Figur 4: Spektrum der Spaltneutronen nach Gleichung (93 1 (ausgezogene Kurve) und Gleichung (95 1 (punkt- 
getrichene Kurve). 
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